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Abstract 

We map adiabatic quantum evolution on the classical Hamiltonian dynamics of a ID gas 
(Pechukas gas) and simulate the latter numerically. This approach turns out to be both insightful 
and numerically efficient, as seen from our example of a CNOT gate simulation. For a general class 
of Hamiltonians we show that the escape probability from the initial state scales no faster than 
|A| 7 , where |A| is the adiabaticity parameter. The scaling exponent for the escape probability is 
7 = 1/2 for all levels, except the edge (bottom and top) ones, where 7 < 1/3. In principle, our 
method can solve arbitrarily large adiabatic quantum Hamiltonians. 
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Standard approaches to quantum computing are based on applying a sequence of unitary 
operations to a multiqubit system, and the solution is encoded in an entangled superposition 
of its eigenstates, which is fragile with respect to decoherence. Regrettably, there is a growing 
realization that this approach is not feasible in the near future. An alternative promising 
paradigm is adiabatic quantum computing 

m 

(AQC), where the solution is encoded 
in the ground state of the system evolving under an adiabatically slow change of a control 
parameter A. Remaining in the ground state protects the system against relaxation and 
dephasing [4|]. 

In principle, any standard quantum circuit can be realized by an AQC . Following; 
the AQC approach, we reach the ground state of a complex Hamiltonian H , which encodes 
the solution to a given quantum algorithm, by adding to H a large bias term, ZH^,, Z ^> 1, 

H(X(t)) = H + X(t)ZH bjJ (1) 

such that H(X = 1) has a non-degenerate, easily achievable ground state. The intrinsic 
imitation imposed on AQC is the finite probability of excitation, via Landau-Zener tunneling 
3,0], at any finite evolution speed (the adiabaticity parameter |A|). This tunneling remains 
even after the effects of external and thermal noise are eliminated. Research on AQC has so 
far concentrated on evaluating and minimi zing the probability of leaving the ground state. 



Both polynomial 6|, |9|, llOj and exponential 111, ll2| slowdown was predicted using AQC, and 



this important question remains open. 

In this paper we investigate what is the probability of escape from the ground state, 
and how far, on average, a system can deviate from the ground state during the adiabatic 
evolution. We do this by mapping the parametric evolution of the system (CQ) on the classical 
Hamiltonian dynamics of a ID gas model with long-range repulsion (Pechukas gas {13I. 
and simulating the latter numerically. This approach, which had not been used so far in 
the field of AQC, turns out to be both physically insightful and numerically efficient way 
to solve dl]). For a general class of Hamiltonians, the probability to stay in the same state 
shows a universal power-law dependence for all the energy levels except for the ground and 
top excited states, and this difference can be qualitatively understood within the Pechukas 
gas model. We also develop a kinetic theory that could, in principle, solve arbitrarily large 
adiabatic Hamiltonians. 

Pechukas gas dynamics. — In this approach we consider the instantaneous eigenstates, 
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\n)(X), and eigenvalues, E n (X), of the Hamiltonian ([T]): H(X)\n) = E n (X)\n). The Hamilto- 
nian H (A) is fully determined by the set of all its instantaneous matrix elements, (m\H (X)\n). 



The latter can be determined from the closed set of differential equations (see Refs. 13l. Il4|): 



, x %m — V m , V m — 2 vai Imn — J]] Irak Ikn \ , \1 I \1 I (^) 

dX dX ^ n {x m -x n y dX \{x m -x k y {x k -x n Yj 

where x n (X) = E n (X), v n (X) = (n\ZHb\n), and l mn (X) = (E m (X) — E n (X))(m\ZHb\n). Equa- 
tions ([2]) describe the classical Hamiltonian dynamics of a ID gas with repulsion, where A 
plays the role of time, and the nth "particle" has a position x n (A) and velocity v n (X). The 
particle-particle repulsion is determined by the "relative angular momenta" / mn (A). 

The mapping of AQC to the classical dynamics of Eq. (T5]) is exact and applies to any 
Hamiltonian. All the information about the Hamiltonian H is contained in the initial 
values, x m (X = 1), v m (X = 1), and l mn (X = 1). Those are given by the appropriate matrix 
elements of the Hamiltonian ([TJ) at A = 1, H(l) = Ho + ZHb = Z(Hb + Z~ 1 Hq). By choice, 
H(l) has non-degenerate, well separated levels, with an easily reachable ground state. The 
initial conditions for Eq. (T2]) can be obtained perturbatively in Z~ x to any accuracy. 

CNOT gate. — Consider a specific example of an AQC. An arbitrary iV-step, M-qubit 
quantum algorithm can be encoded in H using the "ground state quantum computing" 



approach 15], whereby every qubit is represented by an array of 2 x (N + 1) quantum dots 
sharing a single electron. The state of the mth qubit on the (n + l)st step of the algorithm is 
given by the probability amplitude to find the electron on either the quantum dot (m, n, 0) 
or (m, n, 1). The solution is determined by the ground state probability amplitudes on the 
quantum dots (m, N, 0) and (m, N, 1). 

For a simple universal quantum gate, the CNOT (where N = 1, M = 2), Hn can be 

nn 

written as |6|, 1151 ] 

Hcnot = (coioCii — Co 00 Ci )(CiiCoio — CioCooo) + 
( c oiiCii — CooiClo ff:c )(Cii c oii — a x Cioc 00 i) + C'doCooCnCn + CoiCoiCioCio, ^ 

where c' mn j creates an electron on the corresponding dot, and C^ n = (c' mn0 , cj„ nl ). The 
ground state energy of ([3]) is zero. In order to specify the initial state of the qubits before 
the operation, a small correction must be added. 

The results shown in Fig. [T^i correspond to CNOT 1 00) — > 1 00) . To achieve this, the 
term 5H = e(cl 00 c Q00 + cl 00 c wo ), with e = —0.1, was added to H . Remarkably, even 



though the initial conditions for the Pechukas equations (jSJ) were only calculated to first 
order in Z~ x = 0.1, the results agree with the exact diagonalization to four significant 
figures, indicating a high efficiency of the approach. The same accuracy holds for the CNOT 
acting on all the other basis states. All the degeneracies of the Hamiltonian ([3]) are precisely 
reproduced in our approach. Since any ID gas with repulsive interactions naturally expands, 
the decrease of the bias potential XZHb corresponds to a contraction from H to H$. Such 
evolution can be considered a Loschmidtian time reversal [l6] of the natural expansion of 
the Pechukas gas. Note that while the levels generally repel, certain groups of levels can 
cross, due to the symmetries of the Hamiltonian ([3]). 

Statistics of level occupation. — Unlike the simple example of a CNOT gate presented 
above, when considering a general AQC case, we would benefit from the knowledge of the 
statistical behaviour, for a given class of AQC problems (i.e., Hamiltonians Hq). An ensem- 
ble of Hamiltonians H corresponds to a distribution of initial conditions in the Pechukas 
dynamics in Eq. (j2j), over which an appropriate averaging must be taken. We choose a set 
described by one of the random matrix theory (RMT) Gaussian ensembles [l7]. This, in 
particular, means that the distribution of x„(A = l),v n (X = 1) and l mn {\ = 1) is Gaus- 
sian. Such an assumption about the behaviour of a large collection of qubits with varied 



couplings is reasonable and was used recently in 10|, |llL but must still be taken cum grano 
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show that for the 3- SAT problem, 



salis. For example, the numerical calculations in Ref. 
RMT describes well only the bulk of the spectrum. In the case of, e.g., flux qubits threaded 
by a magnetic field, the natural choice is the Gaussian unitary ensemble (GUE) of general 
Hermitian matrices. (If the system has time-reversal symmetry, one should instead use the 
Gaussian orthogonal ensemble (GOE) of real Hamiltonians 171].) Note that the solutions 
of ([2]) are described by RMT only when A — > 0. For most of the evolution, the perturba- 
tion ZHb will impose its own symmetry. Eventually, we are interested in the probabilities, 
P(m\n), for the system to end in the eigenstate \m) at A = after starting in the eigenstate 
|n) at A = 1. 

In the absence of external and thermal noise, the interlevel transitions are due to the 
Landau-Zener tunneling. At each level anticrossing, where two levels approach to a minimum 
separation A min , they may exchange their current occupations: P\(m\n) < — > P\(m + l|n). 
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The probability of this process 
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p LZ = exp [ ^ J] (4) 

V 4nh\(m\ZH b \m+l)\\\\) 

strongly depends on |A|. We always assume a uniform evolution, |A| = 1/T. 

The results of simulations for the GUE are shown in Figs. fTTl2l Remarkably, the probability 
to remain in the initial state by the end of the process scales as T 1 / 2 , with the same power, 
1/2, for all the levels, except the edge ones (top and bottom eigenenergies) . For the edge 
levels the probability scales as T 1//4 in the case of 50 levels, and as ~ T 1 / 3 for the case of 
150 levels. 

The different scaling exponent of the edge levels, compared to the bulk, cannot be only 
due to the long range cubic repulsion in Eq. (T5]) between the levels. The latter is obviously 
responsible for the difference between the bulk and the margin of the spectrum in the number 
of anticrossings (i.e., collisions in the ID gas), but this dependence is a smooth function of 
the level number. It is the edge (top or bottom) position of the level, that is important. 
Indeed, while the separation of the internal levels after a collision is limited by the subsequent 
collisions with the levels on the opposite side, this does not apply to the ground state and 
the top state. This can be clearly seen in the top right inset of Fig. [2^. One also sees there 
that the number of collisions (i.e., avoided crossings) for these edge levels is on average less 
than one. This stresses the special robustness of the ground and top states, precisely due to 
their edge positions. 

Another remarkable result is the behaviour of the r.m.s. deviation of the system from its 
initial state (Fig. Eb). It also scales no faster than a power of T. However, the exponent 
here is a smooth function of the level number n. This means that the observed behaviour is 
not described by a simple diffusion equation. 

Kinetic theory for large adiabatic quantum Hamiltonians . — The above results indicate 
that the approximate ground state energy of a system described by the GUE can be efficiently 
determined by adiabatic quantum evolution. Indeed, to reach the accuracy e, the adiabaticity 
parameter must scale as a power of e. On the other hand, e.g., the time necessary to find 
the ground state energy of a complex system (e.g., a spin glass) with classical annealing 



algorithms generally depends on 1/e exponentially (see Ref. [18[ and references therein). It 
is therefore tempting to conclude that the adiabatic evolution of a quantum computer could 
provide an exponential speedup, on average, for this problem. 



In order to confirm this conjecture and establish the limits of its validity, we need to 
consider larger systems. Then the brute force approach to solving the set ([2]) becomes 
inefficient. We can instead identically rewrite it as a chain of equations (see e.g. [if]]) 
for the microscopic distribution functions Fi(x,v,n) = J2j${x ~ Xj)5(v — Vj)5(n — rij), 
F 2 {x, v, n; y, u, m; I) = J2j,k^( x ~ x j)K v ~ v j)K n ~ n j)Kv ~ Xk)$( u - v k )5(m - n k )8(l - l jk ), 
Gi(l) = Y^jk^i} ~ ljk), ■ ■ ■ Here we included the occupation numbers of the levels in the 
description. After ensemble averaging, these produce the BBGKY (Bogoliubov-Born-Green- 



Kirkwood-Yvon) chain 



161 ] for the averaged distribution functions, fi = {F x ),gi = (Gi) . . .. 



The first equation of this chain is 
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/i(z,u,fi) = 2 q^Y1 J dldydu 



l/l 



; f 2 (x,v,n;y,u,m;l). 



(5) 



(y - x) 3 

Following the standard kinetic approach |16| . the BBGKY hierarchy can be truncated 
by representing the r.h.s. of (j5j) as a term responsible for the long-range interactions (where 
f 2 can be factorized, f2(x,v,n;y,u,m;l) fi(x,v,n)fi(y,u,m)gi(l)), plus the collision 
integral, 1st- This leads to an approximate equation for j\: 



fi(y,u,m)\ d_ 
/ dv 



fi{x,v,n) 



'st, 



(6) 



(y - x) 3 

where the effective repulsion V — J dl g±(l) \l\ 2 , V denotes the principal value of an integral, 
and the collision integral 

J st = 2r st Yj / du (u - v) p LZ (u - v) [/i(a;,u,m)/ 1 (a;,«,n) -/i(2c,u,n)/i(a;,it,m)]. (7) 

Here T St is a constant; the collision integral ([7]) describes the population exchange of two 
anticrossing levels due to Landau-Zener tunneling (j4j). 

If we are only interested in the behaviour of energy levels, but not in their occupation (i.e. 
considering /i = YJ n /i), the kinetic equation for fi will have the form (JS)), but without the 
collision integral. When the full hierarchy of dynamical equations ([2]) becomes intractable 
due to the large size of the system, one can truncate it to only a few lower states. The 
influence of the higher-energy states can be then taken into account statistically using the 
kinetic equation for f\. Due to the long-range interactions in the system, one expects that 
the average behaviour of the low-lying states can be thus accurately predicted. Thus this 
approach could be extended to simulate any large quantum system. 
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Conclusions. — We model the quantum adiabatic evolution of any Hamiltonian by map- 
ping it to the dynamics of a ID gas. Using a CNOT gate as an example, we found that 
the approach is reliable. In the case of GUE Hamiltonians, we found that the behaviour of 
the system is not described by a simple diffusion equation, because the standard deviation 
from the initial state and the probability to stay in it scale differently with the adiabaticity 
parameter |A|. The probability to stay in the same state shows a universal power-law de- 
pendence for all the levels except for the ground and top excited states, and this difference 
can be qualitatively understood within the Pechukas gas model. For the investigation of 
large systems, where direct simulation is impractical, we propose a kinetic approach based 
on the BBGKY chain for the Pechukas gas. The results of this paper indicate that an adia- 
batic quantum evolution can provide an exponential, on average, speedup compared to the 
classical simulated annealing in finding an approximate ground state energy of a complex 
system. 
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FIG. 1: (Color online) (a) CNOT gate simulation for the operation 1 00) — > 1 00) , described by 
the Hamiltonian (P) with Hq = Hqnot + 5H (see Eq. ([3]) and below). Some levels cross due 
to symmetry, but not the ground state. The small orange diamonds on the vertical axis of the 
main panel show the results of direct diagonalization of Hq and coincide with the results of the 
evolution ([2]) to four significant figures, (b) Evolution of 50 energy levels in Eq. ([2]) as a function 
of the parameter A, for a single GUE realization of the Hamiltonian (see text). The energy scale 
is determined by imposing the normalization (AE) = h\X\, where {AE) = (E n +i — E n ) is the 
average level spacing. The r.m.s. amplitude of the fluctuations in E n and {m\Z H^\n) is O.l(AE'). 
The ground state (red) and the highest excited state (orange) are labeled. It is clearly seen in the 
insets (which enlarge the green boxes) that avoided level crossings appear at small values of A (e.g., 
green vertical arrow in the top inset). In the left inset notice the significant separation between 
the ground state and the first excited state. 
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FIG. 2: (Color online) (a) Probability P(n\n) for the system to remain in the initial state \n), as 
a function of 1/T, during the adiabatic evolution (i.e., decreasing the external field with constant 
speed A = 1/T). The data are averaged over 400 GUE realizations with N = 50. The power-law 
dependence P oc T 7 is clearly seen, with two distinct exponents: (i) 7 = 1/2 for all the bulk 
states and (ii) 7 1/3 for the edge states. For example, for N = 50: 7 = 1/2 for n = 2 (49), 
solid (empty) purple circles; n = 5 (45), solid (empty) green triangles; n = 10 (40), solid (empty) 
upturned magenta triangles; n = 25, black diamonds; but 7 = 1/4 for the ground state (n = 1, 
red squares) and the highest excited state (n = 50, empty blue squares). For N = 150 (left inset): 
7 = 1/2 for n = 2 (purple stars), n = 75 (black pluses); but 7 ~ 1/3 for the ground state (n = 1, 
red crosses). However, the average number of avoided level crossings (right inset) is a smooth 
function of the level number for both iV = 50 (red) and N = 150 (black). As expected in the 
absence of external noise, the probability P(n\n) saturates at 1 as 1/T — ► (when 1/T < 0.01). 
(b) The standard deviation ((n — no) 2 ) of the system from the initial state n = no during the 
adiabatic evolution, as a function of 1/T (N = 50). The scaling is approximately power-law, but 
with the exponent smoothly dependent on the initial state. (Inset) Probability P(n\no) to occupy 
level |n) at the end of the evolution, starting from level tiq = 1, 25, 50 (N = 50). Different curves 
correspond (top to bottom peaks) to 10 3 x |A| = 1, 2.5, 5, 10, 25, 50, 75. 



